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ABSTRACT 

iKijak et al] j201 lal) studied the radio spectram of PSR B 1259— 63 in an unique binary with Be 
star LS 2883 and showed that this pulsar undergoes a spectrum evolution due to orbital mo- 
tion. They proposed a qualitative model which explains this evolution. They considered two 
mechanisms that might influence the observed radio emission: free-free absorption and cy- 
clotron resonance. Using the same database we constructed spectra for chosen observing days 
and obtained different types of spectra. Comparing to current classification of pulsar spec- 
tra, there occurs a suggestion that the appearance of various spectra shapes, different from a 
simple power law which is typical for radio pulsars, is possibly caused by envkonmental con- 
ditions around neutron stars. Therefore, the case of B 1259— 63 can be treated as a key factor 
to explain not only the GPS phenomenon observed for the solitary pulsars with interesting 
environments and also another types of spectra (e.g. with break). 
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1 INTRODUCTION 

The shapes of pulsar spectra are very important in understanding of 
the pulsar radio emission. Typical radio spectrum of a pulsar can be 
described using either a simple power law with a negative spectral 
index of — 1 .8 or two power laws with spectral indices of —0.9 and 
—2.2 and a break frequency on average of 1.5 GHz dMaron et al.l 
I2OOOI) . F or several pulsars there were obse r ved a low frequency 
turnover ( ISiebeJigVallMalofeev et al.ll994 . lLorimer et al.lh99a) 
also mentioned about some pulsars having positive spectral index 
or flat spectrum in the frequency range between 400 MHz and 1600 

MHz^ 

iKiiak et alj j2011hl) provided a definite evidence for a new 
type of pulsar radio spectra. These spectra show the maximum 
flux above 1 GHz, while at higher frequencies the spectra look 
like a typical pulsar spectrum. At lower frequencies (below 1 
GHz ), the observed flux decreases, showing a positive spectral in- 
dex toiak et aljr2011bl) . They called these objects the gigahertz- 
peaked spectra (GPS) pulsars. A frequency at which such a spec- 
trum shows the m aximum flux was called the peak frequency. 
iKiiak et al ] J201ld) also indicated that the GPS pulsars are rela- 
tively young objects and they usually adjoin such interesting envi- 
ronments as HII regions or compact pulsar wind nebulae (PWN). 
Additionally some of them seem to be coincident with the known 
but sometimes unidentified X-ray sources from 3rd EGRET Cata- 
logue or HESS observations. 

Candidates for GPS pulsars were selected by iKiiak & MaronI 

amongst those which showed a decreasing flux density at 
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frequencies below 1 GHz. iKiiak et ^ j2007h used multifrequency 
flux measurements for these candidate pulsars and presented the 
first direct evidence of a high frequency turnover. All of them 
have proved to be relatively young and to have high value of 
dispersion measure (DM around 150 pccm^"^ and ab ove). PSR 
B1259— 63 was also inscribed bv iLorimer et al.l d 1995b in the list 
of pulsars with positive spectral indices. This pulsar is in an unique 
binary with a massive main-sequence Be star. B 1259— 63 has a 
short period of 48 ms and a characteristic age of 330 kyr which 
places it amongst middle-aged pulsars. Its average DM is about 
147 pccm~^ and the corresponding distance is about 2.75 kpc. 
Therefore, B 1259— 63 seems a natural candidate to be classified as 
the GPS pulsar. The companion star LS 2883 is a 10-mag mas- 
sive Be star with a mass of about lOMf^ an d a radius of 6Rq 
jjohnston et al ]|l994l) . lNegueruela et al.l ( I20T l[) estimates the mass 
of LS 2883 of about 30Mq, however the result is subject to large 
uncertainties. Be stars are generally believed to have a hot tenu- 
ous polar wind and a cooler high-density equatorial disc. The PSR 
B1259— 63/LS 2883 emits unpulsed non-theiTnal emission over a 
wide range of frequencies rangin g between radio and 7-rays, and 
its flux varies with orbital phase. iKiiak et 'Z H2011al) studied the 
radio spectrum of B 1259— 63 using the available measurements of 
the pulsed flux obtained during three periastron passages (1997, 
2000 and 2004) and showed that this pulsar undergoes a spectrum 
evolution due to its orbital motion. In addition, demonstrated that 
the peak frequency also depends on the orbital phase and therefore 
it varies with the changes of the pulsar environment. They argued 
that the observed variation of the spectra is caused by a combi- 
nation of two effects: the free-free absorption in the stellar wind 
and the cyclotron resonance in the magnetic field. This field is as- 
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Table 1. The results of our best fit choice to B 1259—63 spectra for each 
observing day, "-" and "+" denote days before and after periastron passage 
respectively, year is the epoch of observation. 



Day 


Year 


Best fit 


Comments 


-60 


1997 


turnover 




-48 


1 QQ7 


unclesr 


no trend 


-46 


2000 


broken 




-40 


1997 


turnover 




-24 


1997 


turnover 




-21 




turnover 




-18 


2000 


ILIllHJ VCl 




-h16 


1997 


ILIllHJ VCl 




+20 


1997 


lUllHJ VCi 




-t-21 


2004 


broken/turnover 


sni3.ll positive spectrsl 








itiHpv hplnw C~W^7 

11I\J.Cj\ UCiVJ W £^ vJX IZj 


+24 


2000 


ILIllHJ VCl 




+27 


1997 


lUllHJ VCi 




+28 


2004 


LfUVVCi l£lW 




+29 


2000 


jJUWCl ItlW 




+30 


1997 


1 1 n p 1 f n 1" 

LllJ^iCuJ. 


lIKj LLCJILl 


+33 


2004 


LHJWCl IdW/ UiUJvCll 


Alllall UlCaJS. 


+34 


2000 


UHJJvCll/ ILllllU VCl 


c m ill r\ocil"iA/p CT^pr*ti"Ql 
Alllall UUalLlVC SUCl-llttl 








1 tlHpY r\pl J^\X7 A. / f-rl-T'7 
lllU-CA UCltJW V J 1 1 /. 


-H36 


1997 


turnover 




-H38 


2004 


broken/unclear 


difficult to classify 


-H40 


2000 


unclear 


no trend 


-H42 


1997 


flat 




+47 


1997 


broken/turnover 


small positive spectral 








index below 2.5 GHz 


+55 


2000 


unclear 


no trend 


+63 


2000 


broken 




+64- 


2004 


turnover 




+70 


2000 


power-law/broken 


small break 


-H83 


2000 


broken 




+94 


2004 


broken 




+U3 


2000 


unclear 


significant eiTor bars 


+150 


2004 


broken 




+m 


2004 


power-law 





sociated with the disc and is infused by the relativistic particles of 
the pulsar wind. The apparent resemblance between the B 1259—63 
spectrum and the GPS sugges ts that the same me chanisms should 
be responsible for both cases dKiiak et alj|201 1 J) . Thus, there ap- 
pears a conclusion that the GPS feature should be caused by some 
external factors rather than by the emission mechanism. 

We collected flux density observation s of B 1259— 6 3 mea- 
sured during three periastron passages (IJohnston et al.l Il999l : 
IConnors et al.l I2002I ; I Johnston et al.l l2005h and constructed radio 
spectra for different observing days. Then we analysed their shapes 
in details taking into account the influence of interstellar scintilla- 
tions and expressed implications of B 1259— 63 case on classifica- 
tion of radio pulsar spectra in general. 



2 B1259 63 RADIO SPECTRUM: DETAILED ANALYSIS 

As it was mentioned previously, the typical pulsar spectrum is usu- 
ally described in terms of a simple-power law 



except from a frequency range below some critical value where 
we presume a cut-off which may be caused by a loss of coherence. 
However, there are also known some more complicated spectra, i.e. 
turnover or broken-type. In case of spectrum with so-called break 
the usual next step is fitting two power laws to the inspected data 



Cil/ 



(2) 



where ui, denotes break frequency. Broken-type spectrum is usu- 
ally defined by having two negative spectral indices and becom- 
ing steeper at higher frequencies. Spectrum with turnover, either 
at low or high frequency, c an be modelled in a way proposed by 
iKuzmin & Losovsld fiooH) 
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and Vp = 10^ denotes peak frequency. iKiiak et al] ( l2011ah used 
the same method to analyse both the shapes of GPS and the spec- 
trum evolution of PSR B 1259— 63 for chosen intervals of orbital 
phases and showed that the pulsar undergoes spectrum evolution 
while it orbits around the companion. It is important to mention 
that they used data with available four flux measurements for a 
given observing day except — 18d (i.e. 18 days prior to the peri- 
astron passage), corresponding to the orbital phase which is just 
before the eclipse, with only three observing points. We created a 
database with all available flux density measurements for this pul- 
sar which a llowed us t o re-ex amine the spectrum evolution pre- 
sented by Kiia k et al] l2011al ). Despite using additional observ- 
ing points which weren't been taken into consideration previously, 
it appears that there wasn't significant differences in shapes of 
B 1259— 63 spectra for chosen orbital phase intervals between our 
results and those from Kijak et al. (2011a). However, our analysis 
revealed new implications on classification of radio pulsar spectra. 

To clarify the spectrum evolution, let us examine the 
B 1259— 63 spectrum observed compiled for each observing day us- 
ing the same database. We choose methods which is a combination 
of those us ed by Maron et al. ( 2000) for radio spectra of about 280 
pulsars and lKiiak et alj hoilj) for GPS and B 1259-63 spectrum 
evolution. For reliable results we analyse only spectra with four 
observing points with additional spectrum for — 18d constructed of 
three flux measurements. 



2.1 Data analysis 

To analyse the spectra in details, we attempted to fit both one/two 
power law(s) to the data with errors. Moreover, a simple parabolic 
function was used instead of this given in equation Eq.JSj. We ac- 
complish the fitting using an implementation of the linear least- 
squares fitting algorithm (available in GNU Scientific Library) to 
logarithms of observational data. To choose the best fit we had to 
analyse the results in a complex way. At first, we rejected the evi- 
dent inaccurate approximations and non-physical results, e.g. break 
at frequency below the lowest or above the highest one. Taking into 
account that there are available maximum four data points we could 
not simply compare the values of to denote the best fit. Other- 
wise, we would classify most of spectra as broken-type as fitting 
two power law functions to 2-elements data set each is statistically 
the most accurate one. Let us denote that resting on an existing 
classification of radio pulsars spectra there are three basic types of 
spectra. 

As it was mentioned, we expect broken-type of spectrum to 
have two negative spectral indices. This in turn leads us to classify 
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Figure 1. The fits to the B1259— 63 spectra for the each orbital phase interval, prior (left panel) and after (light panel) periastron passage. 



Table 2. Best fit and parameters fitted to spectra presented on Fig|2] for 
more details see Eq. l[T], (2) and j3); "-" and "+" denote days before and 
after periastron passage respectively, year is the epoch of observation. 



Day 


Year 


Best fit 


Parameters 






+33 


2004 


power-law 


c = 6.9, a = - 


-0.98 




+83 


2000 


broken 


ci = 5.6, ai = 


= -0.18, 










C2 = 69.7, a2 


= -2.3, 


ut = 3.3 GHz 


+20 


1997 


turnover 


a = -4.7, b = 


5.1, c = 


-0.74, 








Up = 3.5 GHz 






+42 


1997 


flat 


c = 3.8, a = - 


-0.15 





spectra showing positive spectral index below some peak frequency 
as turnover spectra. We also specify power-law spectra with very 
small value (close to zero) of spectral index as a flat one. Some 
of spectra showed two different trends (for more details see TabjTll. 
Any other existing cases would be given recognition as unclear. We 
present the results of best fit for each spectrum in Table [T] 



2.2 Results 

iKiiak et alj l201ld) presented detailed study of PSR B 1259-63 
spectra constructed for few orbital epochs. They also analysed pub- 
lished data of the pulsar flux measured during three periastron 
passages in 1997, 2001 and 2004 but they calculated an average 
flux density using only those data with observing points at least 
four given frequencies for a given observing day. We extended the 
database to include all measurements that have been published so 



far. Fig|T] shows the spectra for various orbital phase range which 
can be used as a reference to analysed spectra for given observing 
days. It seems clear that the shape of the PSR B 1259— 63 spec- 
trum depend on the orbital phase. Therefore, we can acknowledge 
that th e PSR B 1259—63 spectrum undergoes evolution. Kijak et alj 
( 12011a ) suggested that free-free absorption due to the stellar wind 
can be responsible for this effect. Moreover, our analysis confirms 
that the shape of spectra are not full y symmetric with respect to the 
periastron point jKiiaketalfeouj) . One can note that flux values 
on the left panel are generally smaller that those on the right panel 
of Fig[T] This effect can be explained by the fact that the waves 
emitted before the periastron p assage trave l longer distance through 
the stellar wind (see Fig. 4 in lKiiak et alj ^201 lah ). Thus, they are 
affected by free-free absoiption stronger than the waves emitted af- 
ter the eclipse. We can also confirm that a peak (or in some cases a 
break) frequency depend on the orbital phase. This can be cause by 
the cy clotron resonance in magnetic field of the disc jKiiak et al.l 
l2011a[) . 

Comparing the results of fitting presented in Tab[T]one can 
note that shapes of spectra for single days are generally consistent 
with those averaged over orbital phase intervals. There also can be 
seen some kind of domination of GPS feature for orbital epoch be- 
fore eclipse. The same time the general scenario leads spectra after 
eclipse from turnover to other shapes. It is important to mention 
that first sign of more typical spectrum (i.e. broken or power-law) 
appears as early as 20 days after periastron passage. The flux at the 
given frequency apparently changes with orbital phases. When the 
pulsar is close to periastron, the flux generally decreases at all ob- 
served frequencies and the most drastic decrease is observed at the 
lowest frequency. Moreover, we noticed all types of radio pulsar 
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Figure 2. The fits to the B 1259— 63 spectra for chosen days ("-" and "+" denote days before and after periastron passage respectively). Each panel 
shows different type of spectra. 



spectra, including a flat spectrum (see Fig. |2)- Close to the peri- 
astron point the spectra of B 1259— 63 resemble those of the GPS 
pulsars. The spectrum for the orbital epochs further from the pe- 
riastron point are more consistent with typical pulsar spectra (i.e. 
power-law and broken). 



3 DISCUSSION 

3.1 Interstellar scintillations influence on B1259 63 spectra 

Pulsars are generally known to be stable radio sources although 
the pulsar signal is affected in a number of ways as it pro pagates 
through the interstellar medium jLorimer & Kramerll2005l) . Thus, 
the measured flux density varies due to diffractive and refractive 
scintillation effect. Scintillation studies have been extremely suc- 
cessful in determining some properties of the interstellar medium 
and the location of the dominant scattering screen. PSR B 1259— 63 
is a object with relatively high dispersion measure (which means 
that its transition frequency is very high) and at the frequencies at 
which it was observed it should be in the so called strong scin- 
tillation regime. This implies that we definitely have to take into 
consideration both refractive (RISS) and diffractive (DISS) scintil- 
lations when analysing spectra for a given day. We used diffrac- 
tive bandwidth A/diss and timescale AfDiss from scintillation 



observations of the pulsar made far from p eriastron at 4.8 GHz 
and 8.4 GHz jMcClure-Griffiths et al.ll 19981) to estimate values of 
these paramet ers at 1.4 GH z and 2.4 GHz assuming AfDiss oc 
/^■^d~°'^ ( Lori mer & Kra mer 2005 ), where / and d denote fre- 
quency and distance respectively. iMcClure-Griffiths et al.l j 19981) 
observed PSR B 1259— 63 during 46 independent observations be- 
tween 1993 August and 1997 February, excluding data taken in 
1993 December and 1994 January during the periastron passage. 
The observations were typically between 30 and 100 minutes long. 
The standard deviation of the diffractive bandwidths and timescales 
values at 4.8 GHz and 8.4 G Hz throughout the orbit is ~30% 
jMcClure-Griffiths et al.ll998h . We estimated the values of A^diss 
to be ranging from 40s at 1.4 GHz to 360s at 8.4 GHz which 
suggests that diffractive scintillations should not affect the aver- 
age flux measurements (observing sessions was usually 4 hours 
long). Roughly estimated refractive timescales vary from 12 hours 
at 8.4 GHz to more than 20 days at 1.4 GHz. However, for lower 
frequencies the modulation index is relatively small which means 
lower uncertainty estimates when measuring flux. However, we can 
definitely point out at least a couple of high frequency observations 
will be affected by refractive scintillations, e.g. for orbital phase 
from -1-1 13 to -1-186 days (see Tab. [TJ including only three spectra, 
one of which is unclear. This leads to conclusion that flux values 
should be averaged over epochs and/or orbital phase intervals to 
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be more reliable. However, it is important to mention that shapes 
of spectra for single days undergo evolution which is similar in a 
significant way to this observed for spectra averaged orbital phase 
intervals. 

3.2 Implications of B1259— 63 case on classification of radio 
pulsar spectra 

iMaron et al.l ilOOd) concluded that the single power law spectrum 
is the most typical for a sample of about 280 pulsars in a wide 
frequency range. Their studies showed that only 10% of pulsars 
requires two power law spectra. Both numbers of low and high fre- 
quency turnover spectra are also still re latively small dSiebeJ 19731 : 
iMaron et al1l20od iKiiak et alj|201 lah . These results indicate the 
possible existence of the "true" shape of spectrum for each pulsar 
and other types of spectra being a rare exceptions of an unknown 
origin. 

It seems clear that the existing classification of radio pulsar 
spectra is stricte morphological and any direct connection between 
extern al or interna l facto rs and the shape of spectrum cannot be 
drawn. iKiiak et al] ( 1201 la) suggested that the GPS feature should 
be caused by some external factors rather than by the emission 
mechanism. One can note that shapes of spectra for single days 
presented in Tab[T] confirm strong environmental influence on the 
pulsar radio emission. This leads us to suggestion that occurrence 
of different types of radio pulsar spectra owes to the environmental 
conditions around the neutron stars. Even assuming that the appear- 
ance of various non-standard spectra shapes in the general popula- 
tion of pulsars can be caused by peculiar environmental conditions, 
it is impossible to isolate different kinds of elements that influ- 
ence pulsar radio emission and analysed them separately. However, 
we can notice a significant difference between typical broken and 
turnover spectrum. 



4 SUMMARY 

We can finalize that detailed analysis of the PSR B 1259— 63 spectra 
for observing days reveals that the shape of the spectrum depends 
on the orbital phase. Moreover, the observed evolution is widely 
consistent with changes of spectra averaged over orbital phase in- 
tervals. Furthermore, the fits presented in TabJTIshow significant 
differences between spectra before and after eclipse. The results 
leads us to conclusion that the PSR B 1259— 63 spectrum varies 
with the changes of the pulsar environment what can be explained 
by combin ation of two effect s: free-free absorption and cyclotron 
resonance (lKi iaketalJl2011ah . At the same time the detailed anal- 
ysis of the PSR B 1259— 63 spectra for observing days revealed ap- 
pearance of all types of spectra. 

Clearly, in the case of this pulsar, various spectra shapes are 
caused by the same physical phenomena, the differences in the 
shapes appear to be of an quantitative origin, rather than qualitative. 
We see the both the GPS type spectrum, and the broken type spectra 
(of various break "strengths", ranging from the flat-normal spectra 
to almost regular power-law ones), apparently appearing due to the 
same effects, only giving a varying amount of absorption. Extend- 
ing that to the study of the spectra of isolated pulsars, we believe 
that the differences between the various pulsar spectra shapes may 
arise from the amounts of certain effects present (like the absorp- 
tion at low observing frequencies), rather than due to the existence 
of different physical phenomena. This would make the classifica- 
tion of the pulsar spectra in general purely morphological, a shape 



description only, and not a tool to discern between various phenom- 
ena causing the spectra to differ from the canonical power-law type. 
Combining this with the study of the GPS pulsars we can conclude, 
that it is most likely the pulsar environment that causes the spectra 
to differ from the standard power-law. In the extreme cases it is 
causing the appearance of a turnover in the spectrum, but in other- 
wise it is changing the spectra only slightly at lower frequencies, 
which results in the detection of a broken-type spectra. 

Therefore we believe, that the case of B 1259— 63 can be 
treated as a key factor to our understanding of not only the GPS 
phenomenon, but also other types of untypical spectra as well. 
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